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SECTION 1 
INTRODUCTION 


The Perkin-Elmer Corporation Optical Technology Division in Danbury, Connec- 
ticut has performed a six month scientific study contract for the NASA/Goddard 
Space Flight Center, Theoretical Astrophysics Branch, Greenbelt, Maryland. 
(Contract No. NAS 5-23327). The study program consisted of the development 
of a computerized mathematical model and the production of simulated ultra- 
violet (UV) spectrograms on film with the aid of the Perkin-Elmer Line Scan 
Image Generator (LSIG), This model was used as a basis to predict the optical 
performance of the International Ultraviolet Explorer (lUE) Scientific 
Instrument which is scheduled for launch in December 1976. 

The two questions of major concern in the study program were, first, the 
effect of the destructive readout process in the SEC vidicon camera on the 
accuracy of spectroradiometric data in the UV spectrogram (or ’’echellogram", 
since it is produced by crossed echelle and concave spherical gratings), and 
second, the question of the extent of overlap between adjacent echelle orders. 

The results of the study indicate that the answer to the first question is 
that, to the extent that linear modeling accurately represents the imaging 
system, the errors in wavelength determination and radlometry introduced 
by destructive readout are not serious, possibly even negligible. The 
second question has not been completely answered. Order overlap is signifi- 
cant at the bottom of the echelle format where the orders are most closely 
spacee. Further work in the area of image enhancement, particularly desmear - 
ing in the cross-order direction, is suggested. 

Section 2 discusses the development of the math model. Section 3 describes 
the preparation of the input echellogram computer tapes which were linearly 
interpolated from the NASA- furnished data. Section 4 gives the details of 
the simulation procedures, Section 5 discusses the results, and gives 
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recommendations for further work. Section 6 contains a brief summary and 
lists the major conclusions* Selected computer printouts of the simulation 
tapes are included in Appendix A, and a derivation of a series expansion for- 
mula for evaluating a Hankel integral transform is presented in Appendix B. 

Positive film transparencies and photographic prints of the simulated echello- 
grams are furnished in Volume II of this report. 
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SECTION 2 

SYSTEM /VWTHEMATICAL MODEL 

2. 1 MODEL SUMMARY AND BLOCK DIAGRAM 

A simplified functional block diagram of the International Ultraviolet 
Explorer (lUE) Scientific Instrument is shown in Figure 1. (A complete 
description of the Scientific Instrument is given in Reference 1. ) 

In this study program we have concerned out selves with analysis of performance 
degradations of only the imaging portion of the system chain extending from 
the telescope through the SEC vidicon camera_, since effects of noise in the 
transmission link and reconstruction artifacts are deemed to be minor consider- 
ations relative to the questions pertaining to the discharge overlap in the 
SEC readout process and loss of resolution in the UV converter. 

In the following paragraphs a linear analytical model suitable for computer 
simulation of the imaging chain effects on UV echelle spectra is developed. 

The destructive readout process is first discussed, and a computer algorithm 
for generating an ’’effective'' point spread function (PSF) for destructive 
readout is presented in Paragraph 2.2. This is followed by a discussion of 
the derivations of the UV converter and guidance (servo jitter) transfer 
functions in Paragraphs 2.3 and 2.4. The echelle spectrograph is assumed 
to make a negligible contribution to performance degradations, as discussed 
in Paragraph 2.5. Cascaded transfer functions and PSF’s are presented in 
Paragraph 2.6, and finally, a discussion of the implementation of the model 
on the computer is given in Paragraph 2.7. The Perkin-Elmer IBM 370/158 
system was used for all the modeling computations. 

2.2 DESTRUCTIVE READOUT OF SEC VIDICON 

We start with the following basic assumptions concerning the destructive 
readout process : 

a. The digitally addressed read beam current profile is circular 
Gaussian (although other profiles may be substituted in the model). 


2-1 




uv 


SEC 


ENCODER 




1 Ultraviolet i'.xplorer (lUE) 

ns t rumen t Functional Block Diagram 


ER-267 












ER-267 


b. The charge removed from a uniformly charged region of the SEC 
target during a single capacitive discharge read operation is proportional 

to the Integral of the product of the charge density Q times the beam current 
profile over the spatial extent of the beam. In other words, a single read 
operation makes a Gaussian-shaped ’•hole'* in the original uniform charge 
density (see Figure 2-b). 

c. Target lag effects are negligible. 

The linearity of the analytical model rests on the latter two assumptions. 

In the following development, we first expand and generalize the one-dimen- 
sional approach outlined in some preliminary calculations done at NASA. 

Next, we make the important simplifying approximation which neglects the 
previous read operations more than a few pixels (3-4) removed from the 
current beam address. This approximation has the very useful result of 
rendering the readout integral space-invariant, which allows us to use the 
powerful convolution theorem for Fourier transforms. The next step is a 
straightforward extension to two dimensions. Finally, the discrete sum 
approximation to the readout integral is made, which is necessary for com- 
puter calculations. 

2.2.1 One -Dimensional Readout Model 

The following is depicted graphically in Figure 2 for a digitally addressed 
Gaussian beam profile spaced by its full -width -half -maximum (FWHM) width. 

The charge readout at the first beam address is given by 

09 

Mf = J_^ QG(x-Xj^)dx, (2.1) 

where Q is the (uniform) charge density on the target, and G(x-x^^) is the 
beam current profile, appropriately normalized, centered at Xj^ (Figure 2-a), 
The charge density remaining at this point is clearly Q[1 -G(x-Xj^ ) ] ( see Figure 
2-b), which becomes the input to the next readout at x^, namely, 

= I' * Q[l-G(x-x^)]G(x-X2)dx. (2.2) 

J -00 
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c) CHARGE REMAINING APTER SIX 



(CHARGE SIGNAL - 6TH R/0) 


Figure 2 . Destructive Readout of Uniformly Charged Target 
(1-dimensional Model) 
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ThuSj, at the k’th beam addresSj, the charge readout is 

p » k-1 

\ Q ^ [l-G(x-x^ ) ] g(x-Xj^) dx 


i=l 


(2.3) 


2.2,2 Space-Invariant Approximation - Effective PSF 

It should be apparent from the foregoing that^ due to the finite extent of 

the readout beam, the previously read pixels more than a few addresses 

(3-4) removed from the current address make negligible contributions to 

the readout integral, Equation (2.3), For example, consider a pixel spacing 

of 1, 5a, where cr is the standard deviation width of the Gaussian function 
2 2 

G(x) = exp(-x flcs ). The value of G four pixels removed from the current 

“8 

position (x = 0 in this example) is exp(-lS) = 1.5x10 , giving 1-G — 1 to 

. 8 

one part in 10 


Hence, only the 1-G factors corresponding to the nearest two or three 
previously read pixels need be retained in the readout integral. Equation 
(2.3). 

Thus, we are able to make the following approximation to the readout integral. 
Equation (2.3): 

* 00 k-1 

J_„ Q ^ [l-(x-x^ )]g(x-Xj^) dx, (2.4) 

i=k-m 


where we have dropped all pixel contributions more than m beam addresses 
removed from the current (k*th) one. 

Next, we introduce the analytically convenient notion of an "effective PSF" 
for destructive readout by associating the product of (1-G) factors with 
G(x-x^), the PSF of the current readout, rather than with the charge density 
Q. This has the computational advantage of allowing the calculation of the 
destructive readout PSF, or its Fourier transform (a complex optical transfer 
function [OTF]), independently of the input echellogram represented by the 
charge distribution Q. Hence, the effective PSF for destructive readout at 
beam address x^ is defined by (see Figure 2-d): 
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^ V k-1 

\ ^ [l-G(x-x^) jG(x-x^)^ (2.5) 

i-k-m 

"where the braces denote the set of x-x^’s over the range of the products. 

If the readout lattice is of uniform spacing, which is the case for the 
SEC vidicon in the lUE instrument, then the G functions are independent of 
k, and therefore, since m is a fixed constant, the expression (2.5) is 
space- invariant. 

Finally, we rewrite the approximation (2.5) to the readout integral, using 
our new definition of effective PSF, and also generalizing to the case of 
a variable charge density Q(x): 

\ -J- ( 2 . 6 ) 

where the notation is to be understood from the definition of P in Equation 
(2.5). This last integral expression exhibits the form of a convolution, 
which, as we have already noted, is important for the purposes of computation. 

2.2.3 Extension to Two Dimensions 

The extension of Equations (2.5) and (2.6) to two dimensions is obvious : 

Q(x) goes over to Q(x,y), and the (x-x^) go over to (x-x^), (y-y^) giving 

({*‘’‘l}'{yyj}> = «,[l-G(x-x^,y-y^)]G(x-Xj^,y-yjj) (2.7) 

J 

and 

~ J-» J ‘^(=‘'y^\,je( '{y-yj} (2.s) 

for the two-dimensional effective PSF and charge readout expressions 
respectively. Since G(x,y) is circular Gaussian, it may be factored into 
a product of the corresponding one-dimensional functions, i.e., 

G(x, y) = G(x) G(y) 
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This will be convenient later on in the construction of the effective PSF 
computer algorithm. Note in Equation (2.7) that we have not specified the 
ranges of the product indices i,j. These will depend upon the particular 
sequence of readout beam addressing employed. In the case under study, a 
two-dimensional square lattice is digitally scanned in raster fashion with 
pixel spacings specified by NASA to be; 

Ax - Ay = 33ui (2.9) 


at the camera focal plane. Therefore, the beam addresses (x^ y^) 
are given by 


x^ = kAx i k = 0, 1, . , , ,N-1 
i = 0, 1, . . .,N-1 


( 2 . 10 ) 


where N is the dimension of the readout lattice. It seems most convenient 
to indicate the (1-G) factors retained in the product of Equation (2,7) 
diagrammatical ly, which we have done in Figure 3, Here we show the previously 
read pixels in the raster which are within a radius approximately 5Ax (or 5Ay) 
of the current readout address The numbers 1-6 indicate the degree 

of nearest neighbor, where 1 is the current beam position. This is helpful 
in keeping track of the (1-G) factors when writing the computer subroutine. 

The number assignments are somewhat arbitrary, but it does not matter, since 
the ordering of the factors in the product is immaterial. A border five 
pixels wide has been indicated in the diagram around three sides of the 
raster within which the effective PSF is not space-invariant, since (1-G) 
factors are progressively lost as the edges of the raster are approached 
(except for the left-hand side where there are no previously read pixels). 
However, very little information is lost by ignoring the edge regions, 
since they comprise a very small fraction of the total format. Moreover, 
the side edges are also where the blaze angle falloff of the spectrograph 
is greatest, so they are likely to be of little interest. 
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2.2»4 Computer Algorithm for Destructive Readout 


With the aid of the labeling sequence in Figure 3-b^ the (1-G) factors may 
readily be determined and the effective PSF constructed. In the space- 
invariant region^ ^ is independent of the beam address y^). Hence^ 
we take the current readout position to be (0, 0), also setting k = 1 = 0 , 
Then, the G functions for the current readout position and previously read 
neighboring pixels are given by (see Figure 3-b): 


Programming Notation t 


G(x)G(y) 

G(x+Ax)G(y) 
G(x+Ax)G(y+Ay) ^ 
G(x)G(y+Ay) 
G(x-Ax)G(y+Ay) 


GXPOS(l)*GY(l) current readout position 
GXP0S(2)*GY(1)1 

GXP0S(2)GY(2) ( nearest neighbors 
GXP0S(1)*GY(2) ( 

GXNEG(1)*GY(2) J 


G(x+2Ax)G(y) ^ GXPOS (3)*GY(1)~J 

G(x+2Ax) G(y+Ay) GXPOS (3)*GY(2) ( second nearest neighbors 


where the nearest, second nearest, etc. neighbors have been written in 
clockwise sequence, starting at the top. 


From this table, the zeroth, nearest neighbor, second nearest neighbor, etc. 
approximations to the effective PSF may be constructed as follows: 


Zeroth approximation (equivalent to nondestructive readout): 
PSF(l) = C*GXPOS(l)*GY(l) 

Nearest neighbor approximation; 

PSF(2) = PSF(1)*(1 - C*OCPOS(2)*GY(1))*(1 - C*GXP0S (2)*GY(2) ) 

*(1 - C*GXPOS(l)*GY(2))*(l - C*GXNEG(2)*GY(2)) 
Second nearest neighbor approximation: 

PSF(3) = PSF(2)*(1 - CVcGXP0S(3)*GY(1))*(1 . C*GXP0S (3)*GY(2) ) 
* 


— 

The notations GXPOS, GXNEG refer to previously read pixels in the 
X > 0, X < 0 half-planes respectively. GY is not so distinguished since 
all pixels are in the y < 0 half-plane. 
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Note that a parameter C (0 ^ C ^ 1) has been introduced which takes account 
of the possibility of less than 100% target discharge at the center of the 
readout beam. 

The above construction has been carried out to the fifth nearest neighbor 
approximation for the circular Gaussian case^ and programmed as FORTRAN 
subroutine SECPSF. A listing can be found in Figure 4, 

2.2.5 Effective PSF Features 

Subroutine SECPSF was used to make perspective plots of the effective PSF 
on Perkin-Elmer * s Calcomp plotter for two ratios of pixel spacing-to-beam 
sigma width: Ax/c = 1.551^ and 2.35 (FWHM spacing). For comparison, the 

circular Gaussian (nondestructive R/0) PSF was also plotted. These are 
shown in Figures 5 to 7. The x- and y-axes have been labelled in units of 
o. The view is from the negative x, negative y quadrant. The PSF maxima 
relative to unit height for the nondestructive case are indicated in the 
graphs. Corresponding numerical printouts of the nondestructive and de- 
structive readout processing arrays appear in Figures 8 and 9. These arrays 
are calculated from subroutine SECPSF at points on the computational sampling 
lattice which is interpolated 2:1 from the readout lattice. (See Section 
2,7 for a discussion of the requirement for interpolation,) 

The most obvious features to be noticed in the perspective plots are the shift 
of the PSF maximum relative to the center of the beam, and the two nearest 
neighbor holes located on the x- and y-axes. The "ripple", or residual 
charge, between these holes is also visible. The scale of the plots is not 
large enough to see the effects of the second and higher nearest neighbors; 
these can be seen in the processing array printouts. PSF profiles for the 
case Ax/ct = 1.551 at constant x and at constant y are plotted in Figure 10 
from the array data in Figure 8. 

The shift of the maximum and the ripple effects are most evident for the 
Ax/ (7 = 1.551 case, which corresponds to the greatest amount of beam over- 
lap. This value is derived from the NASA- furnished parameter values of 
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SUBROUTINE 5ECPSF(X, Y,DX, DY,S| G,C, KNEAR, PSF) 

DIMENSION PSF(l),GXPOS(6),GXNEG(6),GY(6) 

DXX=DX/SIG 

DYY=DY/SIG 

XX=X/SIG 

YY=Y/S|G 

)MAX=KNEAR+1 

DO 1 I =1^ IMAX 

XPSHF=XX+( I -1)*DXX 

XXPSI^F = XPSHF*XPSHF/2 . 

XNSHF=XX-( l-l)*DXX 
XXNSHF=XNSHF*XNSHF/2. 

YSUF=YY+( !~1)*DYY 
YYSHF=YSHF*YSHF/2. 

IF(XXPSHF.GT.25. )XXPSHF^25, 

IF(XXNSHF.GT.25. )XXNSHF=25. 

1F(YYSHF.GT.25. ) YYSHF=25, 

IFCXXPSHF. LT.-25, )XXPSHF=-25. 
lF(XXNSHF.LT.-25. )XXNSHF=-25, 

IF(YYSHF.LT.-25. )YYSMF=-25. 

GXPOSC I ) -EXP( “XXPSHF) 

GXNEG( 1 )=EXP(-XXNSHF) 

1 GY( I) =EXP(-YYSHF) 

PSF( l)=C*GXPOS(l)*GY( 1) 
tF(KNEAR.LE.O)GO TO 2 

PSF( 2)=PSF(1)*(1 .-C*GXPOS(2)*GY(l) )*(l.-C*GXPOS(2)*GY(2) )*{1.-C*GX 
ZPOS( 1)*GY(2) )*(1 ,-C*GXNEGC2)*GY(2) ) 

IF(KNEAR.Ea.l)GO TO 2 

PSF(3>=PSF( 2)*Cl.-C*GXPOSC 3)*GY(1) )*(1,-C*GXP0S(3)*GY(2) )*(1.-C*GX 
ZPOS(2)*GY(3>) 

Z*(1.-C*GXP0SC1)*GYC3) )*(l,-C*GXNEnC2)*GY(3) )*C1."C*GXNEG(3 )»GY( 2) ) 
IF(KNEAR.EQ.2)GO TO 2 

PSF(4)=PSF(3)*C1.“C*GXP0S(4)*GYC1) )*Cl.-C*GXPOS(4)*GY(2) )*(1,-C*GX 
■ ZPOS(3)*GY(3))*C1.-C*GXPOSC2)*GY(4) ) 

Z*(l,-C*GXPOS(l)*GYC4>)*( l.-C*GXNEG(2)*GY(4) )*(1.-C*GXNEG(3)*GY(3) ) 
Z*(1.-*C*GXNEG(4)*GY(2)) 
fF(KNEAR,EQ.3)GO TO 2 

PSF( 5)=PSF(4)*{l.-C*GXPOS(5)*GY(l) ) * C 1 . -C*GX POS ( 5 ) *GY C 2 ) ) * ( 1 . - C*GX 
ZP0S(4)*GY(3))*( l.-C*GXPOSC3)*GY(4) ) 

Z*(l.-C*GXPOS(2)*GY(5))*f l.-C*GXPOS(l)*GY{5))*Cl.-C*GXHEG(2)*.GY( 
Z5)>*(1.-C*GXNEG(3)*GY(4) ) 

Z*(1.~C*GXNEGC4)*GY(3))*(1.-C*GXNEG(5)*GY(2)) 

IF(KNEAR.EQ,4)GO TO 2 

PSF(6)=PSF( 5)*(l."C*GXPOS(6)*GY(l) )*(l.-C*GXPOS(G)*GY(2) )*(1.“C*GX 
ZPOSC5)*GY(3) )*( l.-C*GXPOSC5)*GY(4) )*(l.-C*GXPOS(4)*GY(4) ) 
Z*(l,-C*GXPOS(4)*GY(5))*(l.-C*GXPOS(3)*GY(5))*(l.-C*GXPOS(2)*GY( 
Z6))*(l.-C*GXPOS(l)*GY(G) ) 

Z*(l.“C*GXNEG(2)*nY(6))*(l.-C*GXNEG(3)*GY( 5) )*(1.-C*GXNEG(4)*GY(5) ) 
Z*(1.-C*GXNEG(4)*GY(4)) 

Z*(l .-C*GXNEG(5)*GY(4))*{1.-C*GXNEG( 5)*GY(3) )*(1.-C*GXMEG(6)*GY(2) ) 

2 RETURN 
END 


14:51 : 16 


Figure 4. Computer Subroutine for Destructive 
R/0 of SEC Vidicon 
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Z-AXIS 


+3. 0 a 


Figure 5, Gaussian Beam Current Distribution 
(PSF For Non-Destructive R/O) 
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PIXEL R/0 SPACING: 

Ax = AY = 1.551cr Z-AXIS 



Figure 6. Effective PSF For SEC Vidicon Destructive R/0 
(100 Percent Discharge) 
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PIXEL R/0 SPACING: 
AX = Ay == 2.35a 


Z-AXIS 



Figure 7, Effective PSF for SEC Vidicon Destructive R/O 
(100 Percent Discharge) 
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SECPSF2D TEST 


UX= 16.500000 0Y= 16.500000 XMIN= -S9.00000C XMAX= 99.000000 YMIN= -99.000000 YMAX= 99.000000 

DISCH FRACTION^ 1.000000 DX$EC= 33.000000 DYSEC= 33.000000 SIG= 21.276993 


to 


DR0PSFtPl6) 


0.0 0.0000 0.0 0.0000 0-0 0-0000 0.0 0.0000 0.0 0.0000 0-0 0.0000 0.0 

0-0000 0-0000 0.0000 0.0000 0.0000 0-0000 0.0000 0.0000 0.0000 0.0000 0.0000 0-0000 0.0000 

0.0 0.0000 0.0 0.0000 0.0 O.OOOl 0.0 0-0001 0.0 0.0000 0.0 0.0000 0.0 

0.0000 0.0000 0.0000 O-OCOl 0.0006 0.0016 0.0C16 0.0015 0.0006 0.0001 0.0000 0.0000 0.0000 

0.0 0.0000 0.0 0-0006 0.0 0.0057 0.0 0.0073 0-0 0.0007 0.0 0.0000 0.0 

0.0000 0.0000 0.0001 0.0016 0.0062 0.0388 0.0796 0-0859 0.0307 0.0082 0.0008 0.0001 0-0000 

0.0 0.0000 0.0 0.0026 0.0 0.1027 0.3966 0.3956 0.1696 0-0382 0.0066 0.0003 0-0000 

0,0000 0.0001 0-0008 0.0086 0.0399 0.2106 0.5118 0.6659 0.1977 0-0662 0.0056 0.0006 0.0000 

0.0000 O.OOOl 0.0016 0-0117 0.0565 0.1672 0.2690 0.2150 0-0B88 0.0198 0.0024 0.0002 0.0000 

0.0000 0.0000 0.0005 0.0060 0.0183 0.0667 0.0656 0.0692 0.0200 0.0065 0.0005 0.0000 0.0000 

0.0000 0.0000 0.0001 0-0005 0.0026 0.0060 0.0081 0.0060 0.0026 0.0005 0.0001 0.0000 0.0000 

0.0000 0.0000 0.0000 0*0000 0.0002 0.0004 0.0005 0-0006 0.0002 0.0000,0.0000 0,0000 0,0000 

0.0000 0-0000 0.0000 0,0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 




NOROPSF(Pll)) 
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Figure 8. Processing Spot for 100% Discharge AK = Ay = 1 551cr 
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Figure 9. Processing Spot for 100% Discharge iiX = AY = 2.35a 
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Figure 10. X - and Y - Profiles of Effective PSF 
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Ax = 33|j. and FWHM beam width of 50|j, (=2.35cr)^ and is the value used for the 
first two destructive readout simulations performed in this study. From 
Figure 10 it can be seen that the shift of the peak is about 0.25 pixel in 
the +x direction and about 0.34 pixel in the +y direction. The y shift (along 
the order direction) would lead to an error of about 0. 019 a In wavelength 
determination_, which would appear to be inconsequential, since it is about 
1/5 the wavelength resolution requirement. 


The relative amount of charge readout from a uniformly charged target in a 

single destructive read compared to a nondestructive read may be calculated 

from the ratio of the integral Equation (2.7) of the effective PSF for 

00 


destructive readout IJ P(|x-x^]^ , {y-y^})dxdy to the integral of the 

GO 

destructive readout PSF JJ G(x,y) dxdy. 


non- 


This ratio can be approximated by the ratio of the sums of the respective 
PSF array elements in Figures 8 and 9. For 907= and 1007o discharge and 
Ax/a = 1.551, the ratios are 0.411 and 0.377 respectively. 

It was belatedly realized that the ratio A^/a = 1. 551 is in error, since 
the FWHM width of 50(j,, which was given to us as a "measured beam width", 
should be interpreted as the effective PSF width in the context of our 
math model, not the readout beam width. (The measurement consisted of 
scanning across a very narrow line target, of the order of 5-10|j. in width, 

which gives the one-dimensional effective PSF. ) The above value of A^/cr 

leads to an effective PSF FWHM width of approximately 41(j,, about 20 7, too 
narrow, as can be seen from the PSF profiles in Figure 10. 

The correct value of beam width was found by computing PSF arrays, varying 

the beam width parameter. The results show that a FWHM readout beam width 
of 70.5)1 yields an effective PSF FWHM width of approximately 50)j.. The 
corresponding PSF profiles are shown in Figure 11. Note that the widths 
are different in the two directions, being wider in the x-direction, which 
can be explained by the fewer previous reads in the negative x half-plane 
than in the negative y half -plane (see Figure 3). 
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Figure 11» X - And Y - Profiles Of Effective PSF 
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The peak shift of 0.4 pixel in both directions is somewhat greater and the 
peak amplitude lower than in the previous case because of the increased 
discharge overlap. This PSF was used to rerun the destructive readout 
simulation at \ = 2000A, 

At the request of NASA^ a simulation was also done for an effective readout 
PSF width of loop., (This actually prompted the reexamination of the proper 
definition of PSF width discussed in the preceding.) The profiles for this 
width (in the x-direction) are illustrated in Figure 12. The width asymmetry, 
peak shift and reduction in peak amplitude are even more pronounced. Most 
of the increased width is in the order direction, which may lead to slightly 

a 

reduced resolution within an order. This would be most noticeable at \ =1216A 
where the UV converter performance is best. The increased spread in the cross 
order direction is not enough to significantly increase the order overlap, or, 
indeed, alter the results of the simulations with the 41|i width. This is 
borne out by the simulation results (see Section 5,). Again, the shift is 
mostly in the cross-order direction which does not affect the accuracy of 
wavelength determination. 

2.3 UV CONVERTER 

2.3.1 Introduction 

The NASA- furnished converter data was provided in the form of a square wave 

O o O 

response (SQR) measurements at four wavelengths: \ = 1216A, 1608A, 2000A and 
2500A. Sine wave response (SWR), or MTF, which is necessary for the perform- 
ance calculations, was not directly given but was contained in the "Detector 
MTF" data of Data Sheet 1 of Reference 2, which consists of the cascaded con- 
verter and SEC vidicon MTF’s. Since the latter had been computed assuming 
nondestructive readout of the SEC target, it was necessary to redo the 
detector MTF calculations. 

The calculation of the converter MTF was accomplished in several steps; 

a. Conversion of SQR to SWR via the Coltman formula (Reference 3). 

b. Fitting the resulting SWR curves to an empirical formula 
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Figure 12. X - And Y - Profiles Of Effective PSF 
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M(s) = exp C"(s/s^)^] (2,11) 

which has been shown by Johnson (Reference 4) to represent the MTF’s of a 
wide variety of optical and electro-optical imaging devices by suitable 
choice of the "FTTF index" n and "frequency constant" s^, and 

c. Re-inverting the best fitting "Johnson- formula" MTF curves back 
to SQR as a comparison check with the original SQR data. 

2.3.2 MTF Derivation from Square Wave Response Date 

The NASA- furnished SQR curves for the UV converter at the four wavelengths 
of interest are shown in Figure 13. The dotted portions are extrapolations 
for the purpose of calculating SWR's via the Coltman formula^ which is given 
by 

M(s) = f [ N(s) + N(3s)/3 - N(5s)/5 + N(7s)/7 + N(lls)/ll 

n 

- N(13s)/13 - N(15s)/15 - N(17s)/17 + N (19s)/19“* ‘ ’ J (2.12) 

where M is the SWR at spatial frequency s^ and N is the SQR. At low spatial 
frequencies, the higher square wave harmonics contribute more than at high 
frequencies - hence the need for extrapolation of the SQR data. The SWR's 
thus calculated are shown in Figure 14. (These are the "best fitting" 
curves to the Johnson formula; the corresponding modulation indices and 
frequency constants are also noted in the figure. ) 

2.3.3 Best Fit to Johnson Formula 

The calculated SWR curves were plotted on log-log paper as -f-nM(s) versus s. 

Equation (2.11) thus plotted appears as a straight line. The intercept 

-J^nFbO yields the frequency constant s^, the 1/e point on the SWR curve, 

while the slope is proportional to the modulation index n. "Best fitting" 

(in the sense of matching by eye) straight lines were superimposed on the 

data and s and n determined in the above manner, 
c 
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Figure 14. Sine Wave Response of Proximity Focussed UV 

Converter (Calculated from Square Wave Response 
Data and Fitted to Johnson Formula) 
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2.3.4 Re-Inyerslon to SQR ~ Comparison with Original Data 

The "best fitting" Johnson formula representations of the SWR data shown 
in Figure 14 were re-inverted to SQR data using the familiar sine wave 
expansion of a square wave: 

/ r 

N(s) = j [m(s) - M(3s)/3 + M(5s)/5 - M(7s)/7 + (2,13) 

and compared with the original SQR data of Figure 13, Agreement was within 
17o over the unextrapolated frequency range. 

2.4 GUIDANCE PSF 

2.4.1 NASA-Fumished Specification of Truncated Gaussian Servo Jitter PSF 

The guidance PSF arising from a stellar point source smeared by servo 
jitter has been specified by NASA to be circular Gaussian with q = 0. 5 sec. 

This is truncated at the entrance aperture of the spectrograph. The circular 
aperture of radius = 1.3 sec. was chosen for study. 

At the camera focal plane, the above parameters become o’ = 70. 7u, and r =212. l(j,. 

o 

2.4.2 MTF Derivation 

Boggess (Reference 2) calculated the one-dimensional line spread function in 
the order direction for the truncated Gaussian PSF, which is an error function, 
and also the corresponding MTF. However, the image at the SEC faceplate is 
smeared by the two-dimensional PSF; hence, we redid the calculation. 

It turned out that our result differed only slightly from those of Boggess- 
less than 1% - in the spatial frequency range 0-14 cy/ram (the system MTF is 
dominated by the converter MTF, and drops below 0.01 above 14 cy/mm). Indeed, 
both calculations differ from the Gaussian MTF (no truncation) by less than 
1% in this frequency range. Above approximately 40 cy/mm, the truncated 
Gaussian MTF begins to depart significantly from the Gaussian function, 
reaching zero at 51.75 cy/mm, with "side lobes" above this frequency. 

(Remember that the MTF is equivalent to the Fraunhofer diffraction pattern 
modulus.) Hence, our derivation by Hankel transform is of academic interest 
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only. The work done is briefly summarized below. Additional material is 
given in Appendix A. 


Since the guidance PSF is radially symmetric^ the MTF may be calculated by 
the finite Hankel transform 


M (s) = — ^ 


J (2rrr s) r exp (-r /2 ct ) 
o o o * 


(2.14) 


where is the Bessel function of the first kind and zero order. 

This integral was evaluated by partial integration (Appendix A) obtaining 
the asymptotic expansion 

M (s) 2 : exp(-2irVs^) - exp(-r ^/2a^) S (2.15a) 

o u .p.. ^ ~ V 


n=o 


(2Trr^s)n , 


which gives explicitly the expected Gaussian MTF in the limit r /a»l. 

o 

An alternative derivation found in the literature (Reference 4) gives the 
result 


M (s) 

g 


9 9 £ / O'/ ^ 

exp(r^ /2a ) -1 n=o ° (2xr^s)n+l 


(2.15b) 


which yields the Fraunhofer diffraction pattern of a clear circular 
aperture in the limit r^/a « 1. The latter asymptotic expansion is ob- 
tained also by partial integration but with the u and dv factors inter- 
changed with respect to those chosen in this author's derivation (see 
Appendix B ). 


Both expressions were programmed for numerical evaluation. However_, due 
to limitations in the accuracy of the IBM SSP subroutine employed to cal- 
culate the in the small argtmient region^ this effort was abandoned, and 
direct numerical integration of Equation (2.14) using Simpson’s rule was 
programmed. The resulting radial MTF for guidance smear is illustrated in 
Figure 15. 
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Figure 15, Point Source Guidance MTF 
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2.5 UV SPECTROGRAPH 

The UV spectrograph consists of crossed echelle and concave spherical 
gratings. In this study, the spectrograph optical performance degradations 
were assumed to be negligible. NASA-generated spot diagrams indicate that 
geometrical aberrations in the spectrograph result in image smear that is 
small compared to smear from the other system components, in particular, the 
guidance and UV converter smear. Our experience at Perkin-Elmer with the 
design and fabrication of instruments employing similar grating spectrographs, 
suggests that the off-axis geometrical aberrations may be more serious than 
indicated by the NASA spot diagrams. However, Perkin-Elmer has not investi- 
gated the particular case of the lUE spectrograph design, since it is beyond 
the scope of this study. 

2 

"Blaze angle falloff" or "ripple", the single slit intensity (sine ) func- 
tion, was also to be included in the math model. However, the NASA- furnished 
unsmeared echellogram data includes this in the order direction already, mak- 
ing it unnecessary to add to the Perkin-Elmer model. However, the NASA data 
does not include the ripple factor for the cross-order direction. The sine 
function squared reaches its half -maximum value at ixsp/x — 0, 463jt, where 
s here is the slit width, and p = sin9-sin with 0^ and 0 being the angle 
of incidence and diffraction angle respectively. Thus, at the half-maximum 

point 0.463\/s. For the spherical grating \/s \/d, where 1/d is 

3 -1 

the grating frequency (3.3802x10 cm for the short wavelength grating and 
3 -1 

2.3314x10 cm for the long wavelength grating). At the short and long 
wavelengths of \ = 140 OA, and 2500A, p ^^2 ^ 0.0219 radians and 0.0271 radians, 
which correspond to half-maximum widths at the camera focal plane of 14. 9ram 
and 18.4nim. 

Hence, the falloff in the cross-order direction is not negligible, and 
should be included in a calibration subsystem for the operating lUE hardware 
(or software). We have not included it in the Perkin-Elmer simulation 
model since it is not of central concern in this study. 

2.6 SYSTEM MTF AND PSF 

The component MTF’s discussed in the preceding paragraphs, with the 
exception of the destructive readout MTF (actually the complex optical 
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transfer function^ or OTF^ since the PSF is asymmetric in this case)^ were 
cascaded together for input to the simulation program. They are illustrated 
in Figure 16 for the four converter wavelengths. 

The destructive readout OTF (or PSF, depending on the method of simulation) 
was used as input in an intermediate simulation in a two-step process con- 
sisting of the smear due to the readout process alone, followed by the 
additional smear of the cascaded converter and guidance components. This is 
discussed more fully in the following subsection (Paragraph 2.7). 

The system PSF (again minus the SEC readout) was obtained by Hankel trans- 
forming the cascaded converter and guidance MTF's. These are shown in 
Figure 17 for the four converter wavelengths. These were calculated for the 
purpose of simulation by the direct convolution method but were not used, 
since the frequency domain method was alternatively used, as discussed in 
Paragraph 2.7. 

2.7 COMPUTER IMPLEMENTATION 

The simulations can be carried out by either direct convolution in the 
spatial domain or by linear filtering via Fast Fourier Transform (FFT) in 
the spatial frequency domain, followed by inverse FFT. (The latter process 
is also referred to as "FFT convolution.") Spatial domain simulation was 
the method initially chosen. However, the cost per simulation proved to be 
greater than anticipated to the point of being prohibitive (^$250 per simula- 
tion). Overlooked in the original decision was the need to interpolate the 
data and also the processing spot to be convolved which approximately doubled 
the array dimensions of each. The input echellogram data points had been 
generated on the SEC readout sampling lattice. At all previously read sample 
points, the processing array amplitudes were zero, for 100% discharge. The 
SEC readout effective PSF (see Equation (2.7)) is non-zero between the sampling 
points, but this information is lost in the digital approximation to the con- 
tinuous convolution integral, unless one samples the PSF between the zeroes. 
This means, since the echellogram and PSF must be defined on the same sam- 
pling lattice, that the echellogram must be interpolated to whatever density 


■ ) 
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desired for the sampling of the PSF. Thus^ for example^ at \ = 1216A the 
processing array size^ requirement (truncating the PSF array at « 2Z of the 
peak value) went from 7x7 to 13x13 (see Figure 17). At \ = 2000A, the re- 
quired size for comparable accuracy becomes 17x17. 

Essentially, the digital computation process introduces a sampling in 
addition to that of the readout. (Again, for convenience, we discuss the 
one “dimensional case. Generalization to two dimensions is straightforward 
as in Section 2.2.3.) The convolution integral Equation (2.6) is replaced by 
a sum, i.e., at the k'th beam address 

M-1 

J ” Q(x) ({x-xj)dx ^ ix Z Q(x^) P^({x^-x/|) 

where the x are the M points on the computational lattice (M = 2N-1, 
m 

where N is the dimension of the readout lattice). 

The two-dimensional form of this discrete convolution was performed using 
program "CONVOL" from Perkin-Elmer's Digital Image Processing Programs 
(DIPP) library for the spatial domain computations giving the result of the 
destructive readout process alone, without the other contributions to smear. 

The latter, namely the guidance and converter smear, were simulated by FFT 
convolution of these with the echellogram smeared by destructive readout. 

In the frequency domain, the Fourier transform of the above integral is 
given by the product of the Fourier transforms of Q(x) and P(x), i.e., 

Q(s) and P(s). The discrete transforms are equal to the values of Q(s) 
and P(s) at the points (Ttv=0, L, ...,M-1) on the frequency sampling lattice 
of spacing As which is related to the spatial sampling lattice by As =: 1/MAx. 
The discrete Fourier transform of Q(x, ), denoted as Q (s^), is 

Q(s ) = As E Q(x^) [exp( — ) ] 

m=o 

and similarly for P(a^). 
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In this study, the system PSF was constructed by cascading the component 
MTF's (Fourier transforms of the PSF’s) and then inverse transforming this 
result to obtain the PSF. Hence, the switch to frequency domain processing 
required only the calculation of the echellogram transform. 

Another program package in the DIPP library called "XFORM" was used for the 
frequency domain processing. It consists of a set of subroutines which reads 
the input "scene" (echellogram in our case), performs the two-dimensional 
discrete Fourier transform via FFT, multiplies this by the processing OTF 
(system OTF), which is complex in this case because of the asymmetrical 
destructive readout PSF, performs the inverse transform of the result and 
outputs to LSIG magnetic tape. Optionally, the echellogram transform may be 
saved on tape or disk in case it is desired to process the same scene with 
more than one OTF, which was the case here. To satisfy the discrete Fourier 
transform requirement of periodicity and continuity of the scene at the edges 
(the points and ^), a raised cosine rolloff to zero intensity was applied 
to the echellogram borders extending nine samples in from each edge. 

The echellograms 256 x 400 were processed in 256 x 256 "chunks," overlapping 
by 112 records (order direction) in the center. The spatial frequency sam- 
pling interval and maximum spatial frequency for an N x N image are given 
by 4 F =s 1 /N*^ and F = (N/2)*4F, where ^ is the spatial sampling interval. 
In the case of N = 256, ^ = 16. 5u,, = 0.23674 cy/mm, and F = 30.30303 

cy/mm. 
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SECTION 3 

PREPARATION OF INPUT ECHELLOGRAM 


The NASA furnished tape used in the subsequent simulations contains three 
files: the first two consisting of slit scan spectra, "crowded^' and 

"uncrowded" absorption spectra in the UV, and the third a simulated "un- 
smeared" echellogram made from the former slit scan spectrum by geometri- 
cally transforming the various orders into the echelle format, but without 
smear due to lUE system degradations. 

The echellogram tape was in VICAR format, i.e,, each data point was written 
as one byte. Some problems were encountered in the process of writing a 
program to read this unfamiliar format. These were eventually solved and 
a section of the tape was rewritten in LSIG-compatible format (two 6-bit 
bytes, reversed) and the amplitudes were multiplied by a factor of 32 to 
provide sufficient grey level amplitude for the LSIG. A 150 x 200 section 
from the bottom of the format was taken where the orders are most closely 
spaced. This is illustrated in Figure 18 which is an LSIG reproduction of 
the NASA echellogram. 

The echellogram subsection was next linearly interpolated two-to-one in both 
dimensions to obtain sufficient sample density for the digital convolution 
purposes as discussed in Section 2.6, yielding an array almost four times 
as large, 299 points x 397 records. (One would expect dimensions of (2M-1) 

X (2N-1) in an array generated by linear interpolation from an MxN array — 
in this case, 299 x 399. ) However, the interpolation program terminated 
after 397 records due to a system bug, ) 
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RECORD NO. 275 


RECORD NO. 475 


POINT NO. 600 


Figure 18. NASA Echellogram "REDPSF” (750 x 750) Reproduced 
on LSIG^ Showing Subsection Used for Simulations 
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the relevant parameters of wavelength^ spatial sampling intervals for the 
SEC readout and for computation^ PSF and R/O beam widths_, spatial frequency 
sampling interval^ etc. 

Negative tranparencies of all simulations were made on the LSIG using a 
clear square writing mask. 

The simulation outputs are shown in the deliverable photographs and posi- 
tive transparencies which accompany the report in a separate package. Tape 
dumps of selected subsections of the simulated echellograms are given in 
Figures A-1 to A-7 of Appendix A. 
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SECTION 5 

DISCUSSION OF RESULTS 

5.1 EFFECTS OF DESTRUCTIVE READOUT 

The major results to be noted here are, first, that the wavelength error 
introduced by the pixel shift due to PSF asymmetry is small for the case 
of the 50 micron width. As discussed in Paragraph 2.2.5, the shift of the 
peak of the effective PSF is less than 1/2 pixel in both the order and cross- 
order directions. However, this does not give the whole story. In Figure 
19 a portion of a single order, smeared by the SEC readout alone (Simulation 
2) and integrated across the width of the order (simulating a slit scan) 
taken near the bottom of the echelle format, was plotted in the order direc- 
tion, For comparison, the corresponding portion of the unsmeared input 
(Simulation 1) was also plotted. The abscissa units are pixel numbers for 
the interpolated echellogram (remember that 2 interpolated pixels = 1 readout 
pixel) and the ordinate is grey level (x32 relative to the NASA- furnished 
data). It is evident from this comparison plot that the spectra are nearly 
identical. The SEC readout spectrum is shifted by approximately one read- 
out pixel with respect to the input, whereas the shift of the PSF peak was 
less than 1/2 pixels as noted above. In addition, there is some smoothing 
of the data. The peaks and valleys agree in amplitude to better than one 
percent. Indeed, it appears that If the shift were taken out, the two spectra 
would agree within one percent throughout. We remind the reader that the 
normalized OTF for SEC readout was used for the simulations, and that a gain 
reduction factor of 0.411 (see Paragraph 2.2,5) should be applied to the curve 
for SEC readout. However, this is not a loss of modulation, but merely gain. 

Hence, the only significant result of the destructive readout process for 
spectroradiometric data analysis as predicted by our linear model (assuming 
sufficient noise-free gain elsewhere in the system chain) is the shift of one 
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readout pixel. The smear resulting from the SEC readout does not appreciably 
affect the resolution along the order direction, if the energy spread in the 
cross-order direction is integrated by slit scanning. For the case of the 
100 microns PSF width, the smear has increased enough in the cross-order 
direction to cause a small degree of order overlap in simulation of the SEC 
readout alone (see the printout in Figure A-3). Because of this, and because 
of insufficient time remaining in the contract, the slit scan calculation 
done to obtain Figure 19 was not repeated. However, in this case, one would 
expect a result similar to the 50 microns case, since the overlap is not 
serious. Some loss in modulation should be expected since the PSF width in 
the order direction is doubled. 

5.2 EFFECTS OF GUIDANCE AND UV CONVERTER SMEAR 

These components, when convolved with the echellogram data smeared by the 
SEC vidicon, give rise to the greater part of the total smear, with resulting 
order overlap, which can be seen in both the photographs and the computer 
printouts of Appendix A. We have also plotted representative profiles in 
the cross-order direction near the bottom center of the format. These are 
shown in Figure 20. 

The wavelength dependence can be observed by comparing simulations 8 and 9 
(reference Table 4-1). As one would expect, the greatest smear occurs for 

O 

the \ = 2000A case, where the converter performance is worst. 

The result of doubling the SEC readout PSF width from 50 to 100 microns is 
largely masked by the converter and guidance smear contributions, as can be 
seen by comparing Simulations 4 and 5. 

A peak shift of approximately one readout pixel in the cross -order direction 
can also be observed in Figure 20. 

Order overlap in the complete system simulations is significant enough that 
one would expect a slit scan calculation similar to that done for the SEC 
readout alone (Figure 19) to be in error, particularly for a case integrating 
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Figure 20. Simulation Profiles in Cross-Order Direction at 
Bottom Center of Echelle Format 
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a weak profile adjacent to a strong one. A suggested approach would be to 
try enhancement filtering in the cross-order direction, wherever the overlap 
is serious. 
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SECTION 6 

SUMMARY AND CONCLUSIONS 


This study program has resulted in the generation of a linear mathematical 
model for computer simulation of the optical performance of the lUE Scientific- 
Instrument imaging chain. A series of computer programs has been written to 
implement the model^ and a number of simulated echellograms were generated 
demonstrating the major effects of performance degradations due to the several 
imaging components. In particular^ an algorithm was developed to represent 
the two-dimensional destructive readout process in the SEC vidicon camera by 
an effective PS F which, to a good approximation, is space -invariant. 

The major conclusions of this program are; 

a. Wavelength determination error due to the pixel shift arising from 
destructive readout is predicted to be approximately one readout pixel (33 
micron width), or approximately 0. 06A in wavelength, for the case where the 
effective PSF FWHM width is taken to be 50 microns. , 

b. Order overlap is principally due to the smear of the guidance and 
converter components, and may be a significant problem at the bottom of the 
echelle format, particularly for the case of a weak order profile adjacent 

to a strong one. Enhancement filtering in the cross-order direction may pro- 
vide help in this instance. 

c. Radiometric error due to the discharge overlap of the destructive 
readout is predicted to be negligible, provided that the assumption of a 
linear readout process is a good approximation to reality, and also that there 
is sufficient noise-free gain In the system chain to overcome the loss in 
gain due to destructive readout. Radiometric error may result from aberra- 
tions in the spectrograph (which were ignored in this study) and also in the 
electron beam focus of the SEC vidicon. 
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APPENDIX A 

COMPUTER PRINTOUTS OF SIMULATIONS 1-7 


Figures A-1 through A- 7 are dumps of selected portions of the echellogram 
simulation tapes for simulations 1 through 7 (of Table 4-1) near the bottom 
center of the echelle format. 

The point and record numbers are those of the computational lattice which 
has dimensions of 256 points x 400 records. The lattice spacing in both 
directions is 16,5 microns. The point numbers increase horizontally (cross- 
order direction) from left to rights while the record numbers increase ver- 
tically (order direction) from top to bottom. 

The simulated slit scan discussed in Paragraph 5.1 and shown in Figure 19 
was hand calculated from swaths of data in Figures A-1 and A-2 between point 
numbers 224-241^ records 126-275, 
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Figure A-1. Echellogram Subsection - Tape No, 812 
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SECTION 4 
SIMULATIONS 


The interpolated echellogram made in Task II was first convolved with the 
effective PSF for destructive readout yielding an intermediate echellogram 
on LSIG magnetic tape. At this point, it was decided to switch to FFT 
convolution processing to save money on computer time, which, as we have 
noted in Section 2, was prohibitively expensive for direct convolution 
because of the large size of the PSF processing array required. 

o 

Destructive readout simulations at converter wavelengths of A. = 1216A and 

o 

2000A were completed via FFT convolution using the cascaded MTF’s for the 
converter, and guidance components as inputs. 

This was followed by a nondestructive readout simulation at \ 2000A which 

was done entirely by FFT convolution. 

At a meeting with NASA, the above results were presented^ and plans for 
modifying and limiting the schedule of further simulations were made. It 
was requested by N^A that we do a destructive readout simulation using 
an effective PSF width of 100 microns (double the width used previously) 
because the order overlap appeared from the photographs to be less than 
expected. In the course of selecting the appropriate beam width parameter, 
it was realized that the effective PSF width and the readout beam width 
are not the same, as already discussed in Paragraph 2.2.5. After determin- 
ing the correct beam width parameter for the 100 microns PSF width, it was 
decided to redo the 50 micron case. This was accomplished by FFT convo- 
lution which required the calculation the complex OTF for the SEC readout, 
due to the asymmetry of the PSF. 

A table of simulation data appears in Table 4-1, including the intermediate 
simulations of SEC readout alone and giving titles, figure numbers, and 
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TABLE 4-1 


SIMULATION data 
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Tape 
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Fig. No. 

(Computer 
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ECHFlLQGftftK SUBSECTION 

tape i,2Z - DESTRUCTIVE R/Q S I MUL AT I ON I COMPLE T€ SYSTEM!, LAMROA = 2000 ANGST, SEC PSF WIDTH = 50 MICRONS 


PTS: 

201 

TO 

22 5 

TOTAL 

PIXELS 

* 102400 SUM = 0.1098E 08 
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TO 
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0.10T2E 
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HCHFLLOGRA'< SUBSECTION 


TAPE M2 - DESTRUCTIVE R/U S I NUL AT 1 ON IC OMPLE T E SYSTEM!, LAMBDA = 2000 ANGST, SEC PSE WIDTH = 50 MICRONS 

PTS: 226 TO 250 TOTAL PIXELS = 102400 SUM = 0.1C98E 08 

RFCS: 125 TO 275 AVG = 0.1072E 03 
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FCHELLOGRAf* SUBSECTION 

T4PE HOT - DESTRUCTIVE R/0 S IHUiAT I ONI COMPLE TE SYSTEH), LANB04 = 1216 ANGST, SEC PSF WIDTH = 61 MICRONS 


PTS: 226 TO 250 TOTAL PIXELS * 100866 SUM = 0.1093E 08 

RECS; 125 TO 275 AVG * 0.1086E 03 


45 47 55 73 109 170 251 317 326 272 191 126 81 57 65 60 38 60 67 63 95 166 207 265 230 

44 46 53 70 102 159 233 296 302 252 178 116 76 55 66 39 38 60 69 66 101 156 222 263 266 

64 65 51 66 96 168 216 272 280 236 166 108 72 52 63 30 38 61 50 69 106 166 237 281 266 

63 66 69 63 91 139 202 253 260 217 156 101 68 50 61 38 38 62 51 72 Z12 176 251 297 279 

42 42 48 t>o 86 130 189 236 262 202 166 95 65 68 60 37 38 62 53 75 116 183 261 310 290 

41 61 66 58 8 2 12 3 178 222 22 7 190 136 89 61 66 39 37 38 42 53 76 119 180 269 319 298 

60 60 65 56 78 118 169 211 216 180 128 85 59 65 38 36 37 62 56 77 121 192 276 326 303 

39 39 66 56 76 113 163 202 207 172 123 82 57 66 38 36 37 63 56 78 122 193 275 326 305 

38 39 63 53 76 liO 158 196 201 167 119 80 56 63 37 36 37 62 56 78 122 192 276 325 306 

37 33 42 52 73 108 155 193 197 165 117 78 55 63 37 35 37 62 56 78 121 191 272 322 301 

3b 37 61 52 72 108 156 192 197 166 117 78 56 62 37 35 37 62 56 77 120 189 269 317 297 

3b 37 61 51 72 108 155 193 197 165 117 78 56 62 37 35 37 62 56 76 119 186 265 312 292 

35 3b 41 51 72 109 156 195 200 166 118 79 55 62 37 35 37 62 53 76 118 186 260 306 286 

34 36 61 52 73 llO 159 198 203 169 12C 00 55 63 37 35 37 62 53 75 117 IBl 256 300 280 

34 36 61 52 76 112 162 203 200 173 123 81 56 63 37 36 37 62 53 75 116 179 251 296 2 76 

36 36 61 53 76 115 167 209 216 170 126 83 57 66 38 36 37 62 53 75 115 177 267 237 266 

33 36 42 56 78 119 172 215 220 183 130 86 59 65 38 36 37 62 53 75 115 175 261 278 258 

33 36 63 55 80 123 178 223 228 190 136 88 60 66 39 37 37 62 53 75 1X3 171 236 268 247 

34 37 44 57 83 128 185 232 237 197 139 92 62 67 60 37 38 62 53 76 III 167 226 257 236 

34 33 45 60 08 136 195 263 249 207 146 95 65 49 41 38 38 42 52 73 109 162 217 245 224 

35 39 67 63 93 163 207 258 263 218 153 100 68 50 62 38 38 62 52 72 107 157 207 232 211 

35 49 49 66 99 153 221 275 280 232 163 106 71 52 43 39 39 43 52 72 106 152 199 220 199 

3ft 61 51 71 106 166 238 295 301 268 176 113 75 55 45 40 40 43 53 72 105 149 191 209 188 

37 43 54 75 114 177 257 319 326 268 187 120 79 56 47 41 41 66 53 73 106 166 185 199 177 

38 45 57 61 iZ3 192 279 347 352 29C 202 129 85 61 48 43 62 65 56 73 106 166 179 191 168 

4C 47 61 87 133 209 306 370 383 315 219 139 90 64 51 44 43 45 55 74 105 143 175 183 160 

41 49 65 93 144 227 331 412 418 344 238 150 97 68 53 46 44 47 56 75 1C6 143 172 177 163 

43 52 60 100 156 247 362 451 457 375 258 162 103 72 56 68 65 68 58 77 108 146 170 172 167 

44 54 73 107 lo9 269 394 491 498 408 200 175 llO 76 58 50 67 69 59 80 111 146 17C 169 142 

46 57 77 114 181 290 427 533 540 442 302 188 117 80 61 52 69 51 62 84 116 151 172 160 139 

49 60 81 121 194 312 460 575 583 475 324 200 125 84 64 54 50 53 65 88 123 158 178 170 138 

50 62 35 128 206 333 492 615 623 508 346 213 132 08 67 56 52 56 68 93 130 167 185 174 139 

61 65 39 134 21 6 353 522 653 662 538 365 224 138 92 t>9 58 54 56 72 99 139 1 70 195 180 142 

53 67 92 141 229 371 55C 688 696 566 384' 235 144 96 72 60 56 61 76 106 149 190 206 187 145 

54 69 96 146 238 388 574 719 727 590 399 244 149 99 74 62 59 63 80 113 160 203 218 196 149 

56 71 98 151 246 401 595 744 752 610 412 251 153 102 76 64 60 66 85 121 172 218 232 205 155 

57 72 101 154 253 412 610 763 771 625 421 256 156 IC4 77 65 62 69 90 129 105 234 247 216 160 

58 73 102 157 258 419 620 775 782 633 427 259 158 105 78 66 64 72 94 138 198 250 263 227 166 

5B 74 103 159 260 423 624 779 785 635 428 260 159 105 79 67 65 74 99 146 211 266 278 237 171 

59 75 104 160 261 422 622 775 780 631 425 258 158 105 79 68 67 77 104 154 223 282 292 247 177 

59 75 K-A 159 259 419 615 763 767 620 417 254 156 1C5 79 69 60 79 108 162 235 296 306 257 182 

60 75 104 156 257 412 603 745 748 604 407 248 153 103 79 69 69 01 112 169 246 311 320 266 187 

6C 74 103 157 253 404 586 722 722 582 393 240 149 101 78 60 69 83 115 176 257 324 333 275 191 

59 74 102 155 248 393 566 b93 691 557 376 231 144 99 77 68 70 84 119 182 267 337 345 284 196 

59 73 lOl 152 242 380 544 661 657 529 357 221 13B 96 75 67 70 85 121 187 276 348 355 291 199 

59 73 99 149 236 367 520 626 621 499 330 209 132 92 73 66 69 85 123 191 282 355 362 295 201 

56 72 93 146 230 354 496 594 584 469 318 198 126 09 71 65 68 35 123 192 285 359 365 297 202 

58 71 96 144 224 342 472 560 547 430 297 186 119 85 69 63 67 34 122 192 204 353 363 295 199 

57 73 95 141 219 329 449 526 511 408 277 174 113 81 66 62 66 82 120 108 279 352 356 233 195 

57 69 94 139 214 317 426 494 475 378 257 163 106 77 64 60 64 30 117 183 271 341 345 279 108 

5b 68 93 137 209 307 406 463 442 350 239 152 100 74 61 50 62 78 113 176 259 326 329 266 180 

56 63 92 136 206 298 388 437 412 325 222 142 95 70 59 55 60 75 100 167 246 308 311 251 170 
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APPENDIX B 

EXPANSION FORMULAE FOR HANKEL TRANSFORM OF 
TRUNCATED GAUSSIAN POINT SPREAD FUNCTION 

B,1 INTRODUCTION 

The calculation of the two-dimensional Fourier transform of an axially 
symmetric point spread function (PSF) may be reduced to the Hankel trans- 
form of the radial PSF, In this Appendix we derive an infinite series by 
partial integration for the case of a circular Gaussian profile truncated 
by a clear circular aperture. An alternative formula, which is obtained by 
interchanging the u and dv variables in the partial integration, is quoted 
from the literature (Reference 5). These formulae yield the well known re- 
sults of the diffraction pattern of a clear circular aperture and the 
Gaussian MTF in the limiting cases of < 1 and > > 1, respectively, 

where r^ is the radius of the truncating aperture, and o is the standard de- 
viation of the Gaussian profile. 

B.2 evaluation OF HANKEL TRANSFORM 

We start with the Hankel transform expression, 

r 

M(s) = (1/a^) J J^(2jtrs) exp (-r^/2a^) rdr, (B-1) 

o 

where is the Bessel function of order zero, and s is spatial 
frequency. Letting z j= 2xrs and 6 = 2jt<7s the integral (B-1) becomes 

2 

M ^ (1/6^) I ° J (z) exp (-z^/ 26^) zdz, (B-2) 

0 o 

o 

where 

z = 2jrr s. 
o o 
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We now integrate by parts, letting u=J^(z) and dv = exp(-z /26 ) zdz, and 

using the relation (l/z)(d/dz) (J (z)/z^) = -J 

n n+ i 

Thus 

M = -J^(z)exp(-z^/2^) (J^(z)/z)exp(-z^/26^) zdz. 


Integrating by parts again gives 

z 2 

M = -J^(z)exp(-z^/26^) I ^ (J^ (z) /z)exp(-z^/26^) 

+ 6^ J 0° zdz. 


z 

o 


0 


Continuing this process and evaluating the upper and lower limits of inte- 
gration with the help of the limit *lim(J (z)/z^) = (l/2^nO we get the 

. ” 
infinite series 

M j= 1 - J (z)exp(-z ^/2a^) - 8^[ l/2-(J (z )/z )exp(-z ^/2a^)] 
o o 1 o o o 

+ 5^ [1/2^2! - (J2(*^)/z^^)exp(-z^^/2o^)] 

- 6® [1/2^3! - (Jj(2^)/2^^)exp(-z^^/2a^)] 


which can be separated into the sum of two infinite series 


M = 1 -5^/2 + (1/2!) (6^/2)^ - (1/3!) (6^/2)^ + ... 


-exp(-z [j (z)-6^J (z )/z +6^J (z )/z ^ 

o o ioo zoo 

(z )/z ^ -3 

+ 0 3 o o . J. 
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2 

The first series is recognizable as the expansion of exp(-6 /2) ; hence 
the expansion formula of the Hankel transform is written, finally, as 

M(s) = exp(-2oi:"a^s ) 

OO 

-exp(-r ^/2a^) S (.2jtr s)/(2jrr s)’^, (B-3) 

o n o o 

n=o 


where and 6 have been replaced by their corresponding parameters. 

In the limit ^ > 1, (B-3) clearly reduces to the MTF of a Gaussian PSF. 

It is evident that the series diverges for any fixed for s large enough. 

Hence, the usefulness of the formula is limited to the region 2jras/(r^/a) <1 


B 3 ALTERNATIVE FORllUlA FOR THE REGION 2Tras/(r /a) >1 

o 

In Reference 5 a similar calculation is performed with the u and dv inter- 
changed in the partial integration. The result in our notation is 


M(s) 


exp(r ^/2a^)-l 
o 


” n+1 

2 (r/a) J ^,(2xr s)/(2xr s)" *•’ (B-4) 

o n+ 1 o o 

n-o 


which reduces in the limit r /<7 

o 


< 1 to 


M(s) ^2 J, ( 27 tr s)/(2jtr s), which is the well-known Fraunhofer diffraction 
loo 

pattern for a clear circular aperture illuminated by an incoherent plane 
wave. It is evident that this series converges only in the region 
2nos/(r /a) >1, so that the two formulae (B-3) and (B-4) are complementary. 


B.4 COMMENT 

In practice, the two expansion formulae may provide an aid in analysis, but 
in numerical computations involving the use of a computer it is probably 
faster and more accurate to perform a direct numerical integration of (B-1) 
using the polynomial approximation to the Bessel function and Simpson's 
rule, which was ultimately done in this study. An attempt to numerically 
evaluate (B-3) on the computer using an IBM Scientific Software Package 


\J 
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subroutine for the J ran into difficulty since the subroutine apparently 

n 

is not accurate for the small argument region of the J^, i.e., in the region 

2iTr s < < 1 . 
o 
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